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dt Abstract. Continuous interpolates are described for classical dynamical systems defined 

' by discrete time-steps. Functional conjugation methods play a central role in obtaining the 

, interpolations. The interpolates correspond to particle motion in an underlying potential, V. 

£NJ ' Typically, V has no lower bound and can exhibit switchbacks wherein V changes form when 

turning points are encountered by the particle. The Beverton-Holt and Skellam models of 
i ■ population dynamics, and particular cases of the logistic map arc used to illustrate these 

O j. features. 
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^ ■ 1 Introduction 

The main purpose of this paper is to provide a few new results obtained from the application of 
functional conjugation methods - Schroder theory [T7j - to problems in dynamics. Specifically, 
the interdisciplinary connections of the Goldstone potential [IT] to the population dynamics 
model of Beverton and Holt pQ, and of qualitatively similar potentials to the mathematically 
more intricate population model of Skellam [18], as described below in Section [3j are original, 
so far as I am awaraS- For an illuminating discussion of these and other population models, and 
■ an introduction to the biology literature, see |10j . 

Moreover, while the discrete logistic map, x *— > sx(l — x), has been particularly well-studied, 
especially since the work of Feigenbaum [9], my discussion of the potentials that underlie the 
map is based on a new investigation using Mathematica code to construct series solutions to 
several hundred orders, and then to functionally continue those series. The remarkable and 
interesting features of the model's phase-space trajectories are again novel, and - at least for the 
s = 3 case presented in Section 0]- displayed here for the first time, so far as I know. Although 
the potentials for the logistic map are given by convergent series for s / 1, for s = 1 the series 
for the potential is only asymptotic. Some explicit results are summarized in the Appendix for 
this special s = 1 case. 



*This paper is a contribution to the Proceedings of the Workshop "Supersymmetric Quantum Me- 
chanics and Spectral Design" (July 18-30, 2010, Benasque, Spain). The full collection is available at 



http://www.emis.de/journals/SIGMA/SUSYQM2010.html 



This statement is reinforced by the fact that my colleagues in mathematical ecology and biology were unaware 
of any such connections. 
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2 Formalism 

Since the functional approach is unconventional and not part of the standard syllabus for classical 
mechanics, to make the discussion self-contained and accessible to the reader, it is necessary to 
briefly review some earlier work on this subject, especially that previously presented in [3], SI El E| - 
Of course, the number of known results in the theory of functional methods is considerable, and 
the mathematical literature on the subject is vast. For a thorough survey of those results and 
an excellent guide to that literature, see [13]. 

2.1 The problem 

Given a discrete map (unit time step, discrete scale change, etc.), 
x i-> fi(x), 

how do we interpolate (for intermediate times, scales, etc.)? Such an interpolation is needed 
to interpret this dynamical system as a continuous flow (Hamiltonian, renormalization, etc.) 
sampled discretely. 

As we have emphasized in recent papers [H El [6] , there is an elegant method, based on the 
work of Ernst Schroder [17] . that exploits analyticity in x to find an interpolation in t. 

2.2 The method 

Schroder's functional equation is 
stf(x,s) = 

for some parameter (eigenvalue) s, with to be determined. So f± is the functional conjugation 
of the eigen-scale-parameter s, 

fi(x,s) = (a:, *),*), 

where the inverse function vf -1 obeys Poincare's equation [16] , 

* -1 (aa?,s) = /i(* -1 (x,s),s). 
The nth iterate of the functional equation gives 

s n #(x, s) = . . s),s) s), s), 

with fi acting n times, and thus the nth order functional composition is also given by a functional 
conjugation, 

f n (x, s) = /i(. . . hihix, s),s)...,s) = ^(s^ix, s),s). 

In principle, with \l/ in hand for a particular problem, any functional composition of f\ can be 
constructed. 

2.3 The solution 

A continuous interpolation between the integer points is then given by functional conjugation, 
for any t, 

f t (x,s) = ^ 1 (s t ^(x,s),s). 

This can be a well-behaved and single- valued function of x and t provided that > I /_1 (x,s) is 

a well-behaved, single-valued function of x, even though s) might be, and typically is, 
mufti- valued. 

For generic real values of ^ appearing in the argument of then: 
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• When s is real and s > 1, the t — > oo behavior for ft(x,s) is determined by the large 
argument asymptotic behavior of VP . 

• When s is real and < s < 1, the £ — > oo behavior for ft(x,s) is determined by the 
behavior of VP" 1 near the origin. 

• When s has a phase, ft(x, s) will be complex for generic t. 

However, regarding this last point, although it may well have been of some interest to partici- 
pants at the workshop on Super symmetric Quantum Mechanics and Spectral Design, here I will 
consider only real s > 0. 

2.4 An interpretation 

The interpolation can be envisioned as the trajectory of a particle, evolving from initial x = 
x(t)\ t _ Q . Thus 

x(t) = ft{x,s), 

where the particle is moving under the influence of a potential according to Hamiltonian dy- 
namics. The velocity of the particle is then found by differentiating with respect to t, 



dt 



(lns)s t *(x)(^- 1 (s**(x))) / , 



where any dependence of VP on s is implicitly understood. Therefore, the velocity will inherit 
and exhibit any multi-valued-ness possessed by VP(x). 

2.5 The potential 

Up to an additive constant, in suitable mass units, V(x) may be defined by [U[5], 
V{x) = -v 2 (x), 

and its x dependence follows from that of the velocity profile of the interpolation, 

In s , 
d 7T7T. (2-1) 



v(x) - 



dt 



41n*(x) 



The potential V = —v then emerges as a quadratic in vP/vP', 

^) = -dn S f(|^) 2 , (2.2, 
where any dependence of V or VP on s is again implicit. 
2.6 A functional equation for V 

The potential can also be determined directly from the functional equation [3] it inherits from VP . 
That is, 

V(h(x, s), s) = (J^fifa 8)\ V(x, s). (2.3) 

If the discrete map possesses a fixed point, we may attempt to solve this functional equation 
for V by series in x about that fixed point. 
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In general, if the series can be constructed, it will of course have a finite radius of convergence. 
However, often the series result can then be continued to other x by making use of the functional 
equation (a technique very familiar, e.g., for the T and £ functions) so that accurate values for 
the potential can often be obtained for all x, for any s. Such series and continuation techniques 
may also be used in many, if not all cases, to determine VE 1 and to arbitrary precision, in 
principle. 

2.7 An aside 

For the corresponding supersymmetric systems, 

V±(x) = —v 2 (x) ± iv'(x). 

This is a complex potential, even for real s > 0. Despite the emphasis placed on such systems at 
the workshop on Supersymmetric Quantum Mechanics and Spectral Design, Centro de Ciencias 
de Benasque Pedro Pascual, where this paper was presented, I will not discuss this any further 
here. 



3 Elementary examples 

3.1 The Beverton— Holt model 

Consider the rational map, 
2r 4- 1 

*->/!(*) = (3.1) 
x + 2 

for — 1 < x < 1. This is known as the classic Beverton-Holt (BH) discrete-time population 
model [1], to som^l, in which context the population is proportional to 1 + x. But physicists - 
especially particle physicists - know it by another name based on the underlying potential. To 
see what this is, consider 

s^(x) = m( ^+1) W ith s = 1/3. (3.2) 
\ x + 2 J 

This functional equation is easily solved by series about either fixed point, x* = ±1, to give 
series that are quickly summed.! Choosing x* = +1 results in 

* x = — — p = ~ , 3.3 

x + 1 2 — x 

i / / t / x \ (x + l) + (x- l)(l/3)* , A , 

V{x) = -v 2 {x) = -i(ln3) 2 (l -x 2 ) 2 . (3.5) 

We immediately recognize V(x) as the Goldstone potential p3J [T4"fl . albeit inverted as encoun- 
tered in instanton problems [15] . as shown in Fig. [TJ 

While ^ and ^t -1 as given above have poles, the underlying potential in this case is an entire 
function, so it is a trivial step to consider particle motion not just on the original interval, but 



2 The model was also discussed in an earlier paper on population dynamics: § 4.3 of [18] - the necessary change 
of variables being most readily apparent upon comparing equation (42) in [18] with (|3.4|l to follow. 

3 Schr6der's equation, in particular (|3.2|l . is homogeneous in So one is free to multiply *!/ by a constant. This 
has no effect on v(x) and V(x), as evident in (|2.1|l and (|2,2[) . nor does it change the trajectory, so the underlying 
physics is the same. 



4 The trajectory (|3.4|l is the familiar "S-curve" of logistic growth, wikipedia.org/wiki/Logistic_function 
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Figure 1. An inverted Goldstonc potential. 



on the complete real line, where V is unbounded below. (This is one reason for my chosen title. 
Another will arise in the discussion of the logistic map in the next section.) 

To be thorough, we note that Schroder's method of analysis also works for an inverted quartic 
potential with only one real extremum. However, in this case the relevant discrete map and 
continuous flow would also have purely imaginary fixed points. To understand this, consider 
v(x) = 1 + x 2 with corresponding potential 



V{x) 



;i + * 2 ) 2 . 



(3.6) 



Construction of the zero-energy trajectory, and a \& function and it's inverse, follow immediately 
for the system governed by (|3.6p . Thus 



ty(x) = exp(arctanx), * 1 (x) = tan(lnx), 
x(t) = ^ _1 (e*^(x)) = tan(i + arctanx). 

Similarities to the BH case are highlighted by the discrete time-step map, now given by 

1 + x 



( — + arctan x 



1 



x 



to be compared to (|3.ip . The corresponding functional equation is 



1 



(3.7) 
(3.8) 



(3- 



(3.10) 



and is easily checked to be satisfied by ^ in (|3.7p by use of the addition formula for the tangent. 

It is interesting that the discrete map in f|3.9[) has period four: f^u o f^u o f^u ° f n /4 = id, 
i.e. the identity map, so that f-K/i{f-K/i{f-K/A{f-K/A{ x )))) = x - Naively, this would seem to require 
that the eigenvalue in Schroder's equation must be a fourth root of unity, thus either ±1 or ±i, 
and not e 71 / 4 as given in (]3.10p . This warrants some explanation. 

The solution ^(x) = exp(arctanx) was obtained by solving Schroder's equation for continuous 
time, namely, 



e t ^(x) = 9(x(t)), 



(3.11) 



where x = x(t)\t=o an d where the time-dependent particle trajectory to be used in (|3.1ip is 
the elementary, right-moving, zero-energy solution (|3.8p for the inverted quartic potential (|3.6p . 
However, this potential deepens so rapidly that a particle moving according to (|3.8p always 



3 I thank the anonymous referee for raising this issue. 
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reaches x = +00 in a finite time starting from any initial x > 0, namely, ir=>oo — § — arctan 2;. 
Of course, this time depends on the initial position but it always satisfies iz^oo < \ for any 
initial position satisfying < x < 00. 

Now, at this juncture you could adopt the narrow point of view that it makes no physical 
sense to consider any solution for larger times t > tx=>oo- This would preclude taking t = ir as 
required to produce f n u o f^u o f n u o f^u = f n = id. And that would be the end of our story. 

On the other hand, if you were just a bit more flexible in your thinking, you might consider 
the particle to be moving on a circle, albeit of infinite radius, and identify the point x = +00 
with the point x = —00. Then for time t > t x =>oo the particle's zero-energy motion would 
continue, except that the particle would be coming in from x = —00 and approaching x = 
from the left. But, for zero total energy in the inverted quartic potential (|3.6|) . it would again 
take the particle only a finite time to go from x = —00 to x = 0, namely, t-oo^o = 5, and 
then the particle would continue moving to the right so that it would reach its original, initial 
position in an additional time io=>z — arctan x. 

In this second point of view, for the particle to carry out one circumnavigation of the infinite 
circle and return to its original position, the required total time to is finite, due to the presence 
of the negative, unbounded potential. That is to say, 

to = tx=>oo + t-oo=>0 + to=>x = 7T- 

This is independent of the initial x, and it corresponds precisely to the statement that f w u o 
/tt/4 /tt/4 /tt/4 {x) = U 0*0 = x, the identity map. 

Moreover, since the particle's instantaneous transition from x = +00 to x = —00 is, mathe- 
matically, just a change from one branch of the arctan function to the next, for the solution 
under consideration, Schroder's equation is indeed satisfied when t = ir but only in the following 
sense: 

e^nth branch (x) = *n+lst branch (x), 

where ^ nt h branch^) = exp(arctan nth branch^)) = exp(n7r + arctan prin cipal branch^))- 

Perhaps a better way to describe all this is to project a circle of unit diameter onto a tangent 
line. Then a zero-energy trajectory along the x-axis defined by that tangent, in an inverted 
quartic potential as described above, corresponds precisely to constant angular velocity (oj = 1) 
around the circle. This is the entire content of (I3.8p . In that expression the argument 9 = 
t + arctan x is just the angle between the circle's diameter perpendicular to the tangent and 
a line from the point on the circle farthest from the tangent to the point at position x along 
the tangent line. When re-expressed as a function of that angle, the particle's energy functional 
becomes 

E[x] = {dx/dtf - (1 + x 2 ) 2 E[x(9)] = {{dO/dtf - l)/cos 4 #. 

To conclude this subsection, it is appropriate to point out that both the models described 
here are PT symmetric, i.e. invariant under x 1— >■ —x and t 1— > —t. And indeed, in keeping 
with one of the prevalent themes of the workshop on Supersymmetric Quantum Mechanics and 
Spectral Design, the quartic model's potential (I3.6P differs from the BH model's f)3 . 5f) simply by 
a change of time scale and the complexification: 1 — x 2 — > 1 — (ix) 2 . 

3.2 The Skellam model 

The discrete time-stepped map in this case is 



x 1 — y k(l - e x ), 



(3.12) 
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where A; is a real parameter. For population dynamics, only < x < oo is of interest, but the 
model is well-defined for all the reals (and C as well). Fixed points are x = and x*(k) = 
k + LambertW(— ke~ k ). The nontrivial fixed point x* is real and positive for k > 1, but real 
and negative for < k < 1, when the appropriate branch of the Lambert function is used. For 
k = 1 there is only the trivial fixed point. This is evident from graphs of the map, for various k. 




Figure 2. Skellam map plotted for various k (namely, 1/2 in red, 1 in green, 3/2 in orange, and 2 in 
blue). Dots are the fixed points, x*(fc). 

This model of population dynamics first appeared in equation (39), § 4 of [18]. (After rescaling 
variables: Skellam's THx = x here, and his TH = k. In fact, in his subsequent discussion of 
§ 4.3, Skellam anticipated the Beverton-Holt model by approximating the exponential in (|3.12p 
for small x: exp(— x) « (1 — x/2)/(l + x/2), the difference being j^x 3 — ^x 4 + T^gjX 5 + 0(x 6 ).) 
The various functional equations associated with (|3.12p are: 

k^(x,k) = V(k(l-e- x ),k), (3.13) 
^-\kx, k) = k(l - e-*" 1 ^)), (3.14) 
k 2 e~ 2x V{x,k) = V(k(l-e- x ),k). (3.15) 

It is straightforward to obtain the first few terms in a series solution of these functional equations. 
Unfortunately, unlike the BH case, it is not possible to express the series in closed-form as known 
functions, so far as I am aware. Nevertheless, the series can be used to compute the potentials 
numerically. 

For example, the potential V, expanded about 0, is: 

n*.*) = -(in^-^x3-i (t+ y_ i)2 *' 

+ ifwTilFf ,S + 0|,i) )' (3 ' 16) 

An alternate series solution for the potential, expanded about x*, is: 



+ MFTWTW^ yi + 0(y6) )- (3 ' 17) 



8 



T. Curtright 



where y = x — x* and I = k — x* = — LambertW (— ke~ k ). And yes, the coefficients in these 
two series do have exactly the same form. 

While the functional equation for V does not determine the normalization, the relative nor- 
malization for the two series (|3. 16[) and (|3.17p has been fixed by the requirement that they match 
on the interval < x < x*. With the relative normalization as given, there is also a duality 
between the right-hand sides of (13. 16h and (|3.17j) to the orders shown: x, k <^=^> y, I. Under this 
interchange the expressions are identical in form, as already noted. This equivalence between 
the two series expansions can be established to all orders by inspection of the corresponding 
recursion relations for the series coefficients. 

For all positive k > 1 the potentials for the Skellam model are qualitatively similar to the 
inverted Goldstone potential, just not as symmetrical. Outside the regions between the two real 
fixed points, which are local maxima, the potentials are apparently unbounded below. Between 
the fixed points the potentials have a single local minimum at x m j n (k) ~ x*/2. For any k > 1, 
the leading x 2 and y 2 behavior of V in ()3.16|) and (|3.17p show that an infinite amount of time 
is required for a zero-energy trajectory to reach the lower and upper fixed points, respectively, 
starting from somewhere else. Consequently, the Schroder function constructed about x = 0, 
with \J/(x) = x + O (x 2 ) , will truly diverge to +oo as x — > x*. Another way to say this is that 
will asymptote to a constant as its argument becomes large, lim (z,k) = x*, when k > 1. 

z— >oo 

For positive k < 1, the Skellam model potentials also have structure qualitatively similar 
to the inverted Goldstone potential, except that the nontrivial local maximum and the local 
minimum are at negative x. For k = 1 the potential is best defined as a limit, V(x, 1) = 
limV(x,k) = —jx 4 — yrjx 5 + 0(x 6 ). However, in this latter case the series for V(x, 1) is 

probably not convergent, but asymptotic instead. (We discuss in some detail a similar situation 
for a special case of the logistic map, in the Appendix. The corresponding analysis of the k = 1 
Skellam model is left as an exercise for the reader.) 

We plot the potentials in the region between the two fixed points for a few representative 
values of k > 1, using the series expansions ()3.16j) and ()3.17|) , only extended to tenth order in x 
and y. This order is sufficient for any errors in the graphs to be insignificant, and indiscernible. 




Figure 3. Skellam model potentials, V(x, k), for k = 3/2, 13/8, 7/4, 15/8, and 2 (as successively lower 
curves) from the series about (in green) joined at x*/2 to the series about x* (in orange). Dots on the 
x axis are the fixed points. 

The continuous trajectories for the Skellam model can also be computed numerically, although 
they are not known in simple, closed form. The most direct way to do this is to compute the 
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inverse of the Schroder function from (|3,14p . say as series, and then re- use that functional 
equation to improve upon the numerical results obtained from the series. Since this article is 
mainly about the potentials that underlie the trajectories, I will not present here numerical 
results for 'I' -1 . But I will discuss the series solution of (|3.14p for generic k. 

The functional equation for is interesting in its own right. Let (x, k) = —f (— x, k), 
then (|3,14p becomes 

f(kx,k) = k(e f{x ' k) - 1). (3.18) 
Consider a formal series solution: 

f(x,k)=J2-^- (3-19) 

i=i 

The exponential formula^, 

OO 




B n (bi, . . . , b n ) n 
x , 



n=0 

where B n (b\, . . . , b n ) are the complete multi-variate Bell polynomials^, permits us to write the 
functional equation (|3. 18[) as 

E bjk j x j _ B n (bi,...,b n ) n 

j=l J - n=l n - 

where it is understood that the bj are fc-dependent. Thus the coefficients of x n give the recursion 
relation, 

fc n - 1 6 n = B n (fe 1 ,...,6 n ). (3.20) 

Since all of these complete Bell polynomials depend on their last argument only as B n (bi, . . . , b n ) 
= b n + • • • , with no other dependence on b n , we can rewrite (|3.20p in a technically sweeter form: 

B n (b 1 ,b 2 ,...,b n - 1 ,0). (3.21) 



k n ~ l - 1 



Moreover, since the B n are homogeneous of degree n under the multi-variate rescaling bk — > r k bk, 
without loss of generality we may choose b± = 1, thereby setting the scale of x in the expansion 
(13~T9]M 

Because of the prefactor in (I3.2ip . we see that each b n+ \ is a polynomial in k multiplied 

n 

by a factor 1/ ] [ (k m — 1), and thus singular as k — > 1. More precisely, I find b n +i(k) ~ 
m=l k-+\ 

(fc-i)n ( n + !)! (\Y- So ' writing 
Pn(k) 



b n +i(k) 



n (k m ~ 1) 

m=l 



c http : //en . wikipedia . org/ wiki /Exponent ial_f orm ula] 

'http : / /en . wikipedia . org/wiki/Bell_polynomial#Complete_Bell_polynomials 
8 From (|3.21[) one can show the series (|3.19[) converges for all x when k > 1. 
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the resulting p n {k) are polynomials in k of order n(n — l)/2. With the choice of initial condition 
b\ = 1, these polynomials have only positive integer coefficients. The first few are: 

Pi(k) = l, P2 (k) = 2 + k, p 3 (k) = 6 + 6k + 5k 2 + k 3 , 
p A (k) = 24 + 36k + 46/c 2 + 40/c 3 + 24/c 4 + 9k 5 + /c 6 , 

p 5 (fc) = 120 + 240/c + 390/fc 2 + 480/c 3 + 514£; 4 + 416/c 5 + 301/c 6 + 160/c 7 + 64/c 8 + 14A; 9 + A; 10 , 
p 6 (k) = 720 + 1800A; + 3480/c 2 + 5250/c 3 + 7028/c 4 + 8056/s 5 + 8252/c 6 + 7426fc 7 

+ 5979/c 8 + 4208/c 9 + 2542/c 10 + 1295/c 11 + 504/c 12 + 139/c 13 + 20/c 14 + /c 15 , 

etc. From these, I infer two particular values of the polynomials, for any n, 

T x n 



p n (k)\k=o = n\, Pn{k)\k=i = n!(n + 1)! 

In fact, the coefficients of the lowest two and highest two powers of k are easily determined for 
any re, 

p n (k) =n\ + -(n - l)n\k + . . . + i( n + 2)(n - lJJfenC"- 1 )/ 2 - 1 + fcM"- 1 )/ 2 . 

This result may be established by induction, and probably extended to other coefficients as 
well, given enough time and interest. 

In terms of the p n s, the formal series solution for is 



^ n\ ^ (n + 1) A M , • 

n=l n=l v ' [7 (1 - k % 

i=l 



We may also express this in terms of deformed, fc-bracket factorial^], 

k m - 1 



Nfc ! = II Hfc> [m]k 



k-1 



m=l 

This makes it clear that the series is actually in terms of the rescaled variable x/(l — k). Thus, 



^ (x,k) = x + x y j 



^(re + 1)! [n] k \ 
As a special case, this has a simple limit as k — > 0, namely, 

lim (y-\x, k))=x + xY -. -r Pn{ )X =x + xY x n = - ln(l - x). 

k^o K y " ^Wre + l)! 1 ^re + l v ; 

n=l n=l 

Another special case yields a direct connection to the BH model, also by taking a limit. Again 
interchanging limit with summation we find 

1 p n {l)x n 



( 1 \ 1 



■;?! 



' x \ n 2x 



.2/ 2-x 

n=l 



where ^ B jj is the inverse Schroder function in (|3.3p . 

"http: //en.wikipedia. org/wiki/Q-bracket#Relationship_to_the_q-bracket_and_the_q-binomial 
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A numerical study of 'I' and the corresponding continuous trajectories for the Skellam model 
for general k will have to be presented elsewhere. For k 7^ 1, the trajectories are qualitatively the 
same as the "S-curves" of the BH model, as can be inferred just from the qualitative similarity 
of the inverted Goldstone potential to the Skellam model potentials, such as those exhibited in 
the Figure above. 

Further features of the functional formalism will now be explored for a different, perhaps 
better-known class of models. 

4 The logistic map 

A well-studied and interesting map is given by [21 [7] 
x i — y sx(l — x), 

on the unit interval, x 6 [0, 1], for parameter values s £ [0, 4]. The maximum of the map is s/4, 
obtained from x = 1/2, so without loss of any essential features, for a given s, we need only 
consider x S [0, s/4]. The map has fixed points at x = and at x* = 1 — lis. 
Iteration of this map shows a wealth of interesting, well-known featurea 10 !. 

• With < s < 1, x n — > independent of the initial x. 

n— >oo 

• With 1 < s < 2, x n 1 — -, monotonically. 

• With 2 < s < 3, x n — > 1 , oscillatorially. 

• With 3 < s < 1 + \/6 « 3.45, x n oscillates between two (s dependent) values. 

• With 3.45 < s < 3.54, x n oscillates between four (s dependent) values. Etc. - the so-called 
"period-doubling cascade" . 

• At s « 3.57 the period-doubling cascade ends, and the x n sequence becomes chaotic. Slight 
variations in the initial x yield dramatically different results for the sequence of x n . 

• Chaotic behavior holds for almost all larger s up to and including s = 4. 
Schroder's equation for the logistic map is 

s^(x, s) = ^f(sx(l — x), s). 

The inverse function satisfies the corresponding Poincare equation, 

^/~ 1 (sx,s) = s* _1 (x,s)(l - ^ _1 (x,s)). 

As originally obtained by Schroder, there are three closed-form solutions known, for s = —2, 2, 
and 4 (s = —2 is related to s = 4 by functional conjugation): 

*(x, 2) = - i ln(l - 2x), V-^x, 2) = I (l - e~ 2x ) , 
^(x,4) = ( arcsin -y/x) 2 , ^ _1 (x,4) = (sin-v/x) 2 . 

Note that ^(x,4) is multi-valued, but ^ _1 (x,4) is single-valued, for < x < 1. 

While the solutions for \I/ for other values of s are not known in a simple closed form, they 
may be constructed as convergent series. The functional method then reveals that there are well- 
defined expressions for the continuous time interpolation of the logistic map for all s values. 

1[ 'http: //en.wikipedia. org/wiki/Logistic_map 
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Moreover, the potentials needed to produce the explicit, continuously evolving trajectories of 
the envisioned particle, for all s, can now be obtained. These potentials had not been obtained 
prior to our work, so far as we can tell from the literature. For s > 3 these potentials are 
unbounded below in an interesting way. 

A new feature for the potentials so obtained is that they show the "switchback effect": If 
and when the particle reaches a turning point, the potential changes form. This requires a bit of 
explanation. 



4.1 The s = 4 case in more detail 

The trajectory, velocity, and potential are given in this case by 
x(t) = * _1 (4 t *(x,4),4) = ( sin (2* arcsin x/x~)) 2 , 

= (\a.A)yJx{t){\ - x(t))arcsin ^/xjt), 



V(x(t)) = (ln4) 2 x(t)(x(t) - l)arcsin 2 yfxjt). 

The last expression appears to be time-translationally invariant. However, the velocity function 
has branch points (i.e. turning points) at x(t) = and x(t) = 1, so some care is needed to 
determine which branch of the potential function is actually in effect, at any given time, since 
the turning points are encountered at finite times. 

Initially, the particle is clearly moving in a potential where all the functions are principal 
valued for x > 0, 



Vq(x) = (ln4) 2 x(x — 1) arcsin 2 \fx. 



(4.1) 
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Figure 4. The primary potential, Vq(x). 

This potential is again unbounded below, as x — > — oo, but the structure for x > is more 
interesting. So, consider the particle to be moving to the right, with initial x > 0. 

Since the RHS of Vo does not have a reflection-free, "infinity pool" edge, the zero-energy 
particle initially moving to the right will reach x = 1 in a finite time, 



At (x) 



1 

hoi 



dy 



VW - V) ( arcsin y/y) 



1 

m~2 



In 



tt/2 



arcsin \ x 



It turns around at this time, but its subsequent motion is that of a zero-energy particle moving 
to the left in a modified potential, V\{x). This follows just from inspection of v(t), and setting 
V\(x) = —v 2 (t) after the encounter with the turning point. 
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This modified potential does not have reflection-free edges at either side, as evident in its 
explicit expression, 

Vi(x) = (}n4) 2 x(x — 1)( — 7r + arcsin yfx) 2 . 




Figure 5. The first switchback potential, V\{x). 

Moving in this modified negative potential, it now takes the zero-energy particle a finite 
amount of time to travel from x = 1 down to x = 0, as given by 

Aii = 1. 

Upon reaching x = 0, the particle is again reflected, and the exact solution x(t) exhibits another 
change in velocity profile that keeps E = 0, but again alters the potential for the return trip 
towards x = 1. The potential becomes 

Vi (x) = (ln4) 2 x (x — 1) (it + arcsin ^/x) 2 . 




Figure 6. The second switchback potential, Vi{x). 
Notice that the sequence Vq — > V\ — > V2 involves successively deeper, negative potentials. 
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And so it goes. The motion of the particle between its encounters with the Pth and (P + l)st 
turning points is determined by the potential: 



V P (x) = (ln4) 2 x(x- 1) (-1) 



1 + P 



7r + arcsin W x 



where [• • • J is the floor function and arcsin is understood to take principal values. That is to 
say, the potential deepens as P increases. 

The corresponding velocity profile speeds up: 



v P (x) = (ln4) v / x(l-x) ((-1) J 



1 + P 



-/r + arcsm Wx 



The particle will either be traveling to the left, with vp{x) < for odd P, or traveling to the 
right, with vp(x) > for even P, with its speed increasing with P. This effect is clearly seen 
upon viewing animations of the s = 4 traiectorieJ^l. 

It is instructive to plot E(t) = v 2 (x(t)) + Vp(x(t)) versus t for various P, to check E = 0, as 
well as to see how the energy would not be conserved if the potentials were not switched. 




Figure 7. E(t) for initial x = 1/2, using potentials Vp, with P = 0, 1, 2, 3, and 4, as shown in green, 
orange, red, purple, and blue, resp. 



For particles with initial v > 0, another way to picture the zig-zag motion of the particle is 
in terms of the total distance traveled. In this point of view, the successive Vp patch together 
to form a continuous potential V{X) on the real half-line, X > 0. Indeed, the half-line may 
be thought of as a covering manifold of the unit interval in x, with the previous multi- valued 
functions of x now single- valued functions of X. In this way of looking at the particle's motion, 
the potential is again unbounded below (with a downward parabola as envelope), even though 
at any time, x > 0. 

Alternatively, the motion may be visualized as a trajectory on the sheets of a Riemann surface. 
Consider the particle moving on the complex x plane, and not just the real line segment [0, 1] , 
the endpoints of which are now branchpoints. There are cuts from +1 to +oo and from to — oo. 
The particle first moves along the real axis, approaches the cut at +1, and then goes around 
the branchpoint such that 

\fx — > y/x, \J\ — x — > — y/l — x, arcsin \fx — > ir — arcsin ^/x. 

The particle returns along the real axis to the origin and encircles the branchpoint at 0, such 
that 

\fx — > — \/x, Vl — x — > y/l — x, arcsin \fx — > — arcsin \fx . 

n http: //server .physics .miami . edu/~ curtright/Schroeder .html 
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V(X) 



-10 




1 2 3 X 

Figure 8. V(X) where the total distance traveled by the particle is X — x(t)\t=o- 



The particle then goes back to +1 and goes around it once more such that again 
y/l — x — > —y/l — x, 



arcsm \/x — > tt — arcsinvs. 



The trajectory continues in this way, flipping signs and adding tts according to the formulae 
for Vp and vp. 

4.2 Series solutions for other values of s 

The functional equation for the potential that underlies the logistic map is 

V(sx(l - x), s) = s 2 (l - 2x) 2 V(x, s). 

Applying the method of series solution about x = 0, with initial conditions V(0, s) = 0, 
V'(0, s) = 0, and V"(0, s) = -2 In 2 s, we find: 



V(x,s) = -(In 2 s)U(x,s), 

/ oo 

U(x,s) =x 2 ( 1 + J^a„(s) 



(4.2) 



ai 



n=l 



l-S 



5 - 3s 



( a -l)2( S + l) : 



provided s ^ 1. The higher coefficients in the expansion are determined recursively by 

1 i+a 



1 



(1 - S n + 2 ) 



Aa n+1 -Aa n + (-l) n j aj i 



V 



n + 2 - j / 



for n > 1, where [• • • J is the floor function. This series solves the functional equation for any s, 
within its radius of convergence, R(s). From numerical studies, that radius depends on s as 
follows: 



R(s) 



lim sup (lan(s)l 1 /™) 



1 

2 

1 



s 
k 4 



if < s < -, 
2 

if - < s < 2, 
3 ~ 

if 2 < s < 4. 



(4.3) 



In particular, i?(4) = 1, as is evident from (|4.ip . while taking (|4.3p at face value, -R(l) = 0, 
suggesting an asymptotic series. More on the s = 1 case as a limit of (|4.2p can be found in the 
Appendix. 
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4.3 Functional continuation of the series solutions 

The sequence of switchback potentials, i.e. the various branches of the analytic potential function, 
can be obtained from the functional equation for the potential. This follows from re-expressing 
the functional equation as 

U±(x) = s(s - 4x)U(x±), x± = - ± ^V 1 ~ 4x/s. (4.4) 

One of these potentials (U-) reproduces the original series for the "primary potential" expanded 
about x = 0, as may be seen by direct comparison of the series or by numerical evaluation. 
Alternatively, the other (t/+) gives the potential on another sheet of the function's Riemann 
surface, hence the first "switchback potential" for s > 2. This is easily checked against the 
closed-form results for the s = 4 case. 

Since U- has built into it zeroes at both x = and x max = s/4, it is actually a more useful 
form than just the series about x = 0. When 2 < s < 4, the series has radius of convergence 
R(s) = s/4. However, from using £7_ instead of the direct series, we have convergence over the 
whole closed interval, x £ [0, s/4], since thenO < \-\\J\ -4x/s < 1/2 < R(s) when 2 < s < 4. 
Thus we need only evaluate U appearing in U- within the region of convergence of its series 
about zero. That is to say, first construct the series for U, and then from that series build U-. 
Finally, identify this with Uq, the primary potential in the sequence. 

Similarly, may be identified with the first switchback in the sequence, U\, but for better 
convergence properties it is useful to build {/+ from Uq instead of U. By doing this when 
2 < s < 4, we have convergence over the whole closed sub-interval, x £ [jqS 2 (4 — s) , s/4] , with 
zeroes of U\ built-in at the end-points of the interval. 

This process may be continued indefinitely, through successive application of the basic sub- 
stitution U{x) — > s(s — 4x)U(x±). For example, the next set of potentials in the sequence is 
U±±±(x) = s(s — 4x)U±±(x±), etc. In general, the nth iteration of the procedure gives 

U±± . . . ±(x) = s(s - 4x)U± ±... ±(x±). (4.5) 

n times n — 1 times 

Finally, at each iteration, we must select appropriate switchback potentials out of the 2 n different 
expressions. In particular, we note that many of the U±±...±(x) will be complex- valued for the x 
intervals under consideration, and therefore they are not of immediate interest since they do not 
govern the particle's evolution along the real axis. (The continuation of the particle trajectory 
into the complex plane is outside the scope of this paper.) 



4.4 Numerical examples 

Consider the specific case s = 5/2. This is representative for 2 < s < 3. For this case, the 
turning points converge onto the nontrivial fixed point x* = l — 1/s = 3/5. This is evident in 
the following graphs. For s < 3, the potential sequence is actually bounded below, for positive x 
(but not for negative x). In fact, for s < 3 the potentials in the sequence become progressively 
more shallow for x > and tend to disappear! 

The turning points and the behavior of the potentials near =3/5 are shown more clearly 
below, for the first three potential switchbacks. 

As the zero-energy particle moves through this sequence of increasingly shallow, narrowing 
potentials, its average speed decreases, giving the appearance of a dissipative system. Never- 
theless, even as the particle motion subsides upon convergence into the fixed point at x = 3/5, 
energy is rigorously conserved through changes in the potential. 

For other values of s, the sequence of potentials can also be accurately determined using 
series solutions augmented by functional methods. 
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Figure 10. Logistic map potentials for s = 5/2, near the fixed point = 3/5. 



Rather than plot the potentials, another useful way to envision the dynamics is to construct a 
phase-space trajectory [5]. This is easily done for those cases where the trajectories are known in 
closed-form. For cases where numerical analysis is required, it is also relatively straightforward to 
construct graphs of the phase-space curves using the series solutions, augmented by continuation 
through use of the functional equations, although for the phase space the relevant functional 
equation is that for the momentum (note that p = 2v for the mass units we have chosen). This 
is the square-root of the potential equation (|2.3h . namely, 

p(fl(x,s),s) = f-^fi(x,s)Jp(x,s). 

(In the mathematical literature, an equation of this form is known as a "Julia functional equa- 
tion" \12\ 113].) The resulting momentum is multi- valued, in general, with the various branches 
given by analogues of the recursion relation (|4.5p . appropriately signed. That is to say, 

p±±...±(x) = ±y/ s(s - Ax)p± ± . . . ±(x±), 



n times 



n — 1 times 
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with x± as defined in ()4.4p . Some of these are actually complex- valued, so for real trajectories, 
it is necessary to choose from among the 2 n possibilities at each stage of the recursion. (See the 
discussion in [3], especially for the s = 10/3 case.) For 2 < s < 3, and for s = 4, the complete 
infinite sequence of real-valued momentum branches is given recursively by 

Pn{x) = -i/s(s-4x)p n _i Q + \^J\ - Ax/s\ . 

For 3 < s < 4, it is necessary to consider additional branches from among the p±±—± (x). 

I illustrate this by plotting the first fifteen branches of a zero-energy phase-space trajectory 
for the marginal s = 3 case. In some ways this is representative of all s > 2 cases. Among other 
things, the plot shows that the system is quasi- Hamiltonian [5] since the resulting trajectory is 
not single- valued. On the other hand, the trajectory possesses such elegant features that I cannot 
help but make the suggestion that the underlying analytic potential is somehow fundamental, 
perhaps even expressible known function. 




Figure 11. Phase-space trajectory for s = 3, as given by dx/dt versus x. As t advances, the trajectory 
converges horizontally, but not vertically, towards x* = |. 

One feature that stands out in the graph is the value of p n {x*) for = 2/3. From the 
recursion relation in this s = 3 case, 

p n (x) = - v / 3(3-4x)p n _i Q + ^l-ix/^j , 

it follows that p n (2/3) = — p n _i (2/3), exactly. Thus, while the phase-space trajectory contracts 
horizontally, towards the asymptotic fixed point at x* = 2/3, the vertical height at x* = 2/3 is 
fixed and nonzero. Numerically, we find \dx/dt\ x=2 / 3 = 0.291464. Remarkably, the transit times 
for the corresponding V n >i branches of the potential are all At = 1, even though those branches 
also shrink to zero width as n — > oo, but with fixed depth, V n (2/3) = —0.0849511. 

For additional details, and numerical results, please see [H [5j [3]. However, I do not wish 
to give the impression that there is nothing further to be done here. To determine the large 
time behavior of maps, there is much to be done. While the asymptotic form of can be 
determined for non-chaotic values of s, I do not know what it is for generic, chaotic versions 
of the logistic map. The one chaotic case where it is known exactly, in simple closed- form 
(i.e. s = 4), is easily understood, but raises many questions about the behavior for other nearby 
values of s. 
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A Series for the s = 1 logistic map potential 

For this case the appropriate limit of (|4.2p gives 



25 



V(x, 1) = lim V(x, s) = -x* - 2x b - Ax" - —x 
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and solve by iteration the functional equation that should be obeyed by V(x, 1), namely, 

V(x(l - x), 1) = (1 - 2x) 2 V(x, 1). 

If the formal series (jA.ip is constructed to O (x 25 ) or so, it appears to have a radius of convergence 
of R ~ 1/2. But, this convergence is seen to be illusory to higher orders. Different behavior sets 
in around O (x 30 ), where the successive |c n | 1//n used in the lim sup determination of R, (|4.3p . 
begin to grow linearly, on average, for n > 30, as shown here (purple curve, with wiggles). 
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Figure 12. Icnl 1 /" for s = 1 and 1 < n < 500. The blue horizontal is 2, while the sienna line is 

( 2 -«/2 e -3n/2 n! )V" 

That is to say, \c n \ ~ L n exp (n In n) for large n, with L w w ^ 5/2 = 5.804 29 x 10 -2 . 

This immediately brings to mind Stirling's formula, n\ ~ \/2vrn (-)", and in fact, a comparison 
of c n with /„ = 2~ n / 2 e~ 3n / 2 n! is striking, as shown in Fig. [12j (The most important feature in 
the figure is agreement between the averaged slopes. A better overall fit to the c n , on average, 
is achieved by shifting the sienna line slightly to the left, for example by including an additional 
yfn factor in f n .) 

All this is compelling evidence that (|A.1|) itself is not convergent, but rather an asymptotic 
series. This is consistent with general results in the mathematical literature [H[T3]. Indeed, 
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the f n "fit" to c n can be obtained from an asymptotic approximation of a simple integral: 
s f°° v 1 ^2 f 3x - 2V2e§ \ . / 1 r 3 \ 

= / e ~ V E ^ n 2- n/2 e- 3n/2 2/^y ~ ^ 

n=0 n=0 

Taking the Cauchy principal value for X{x) gives finite numerical results for all x. These results 
might be useful to compute corrections to the polynomial approximation for V(x, 1) as con- 
structed from truncating the series (jA.ip . I leave further analysis of this interesting but peculiar 
case to the reader. 
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